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Entanglement plays a central role in our understanding of quantum many body physics, and is 
fundamental in characterising quantum phases and quantum phase transitions. Developing proto¬ 
cols to detect and quantify entanglement of many-particle quantum states is thus a key challenge for 
present experiments. Here, we show that the quantum Fisher information, representing a witness for 
genuinely multipartite entanglement, becomes measurable for thermal ensembles via the dynamic 
susceptibility, i.e., with resources readily available in present cold atomic gas and condensed-matter 
experiments. This moreover establishes a fundamental connection between multipartite entangle¬ 
ment and many-body correlations contained in response functions, with profound implications close 
to quantum phase transitions. There, the quantum Fisher information becomes universal, allowing 
us to identify strongly entangled phase transitions with a divergent multipartiteness of entanglement. 

We illustrate our framework using paradigmatic quantum Ising models, and point out potential sig¬ 
natures in optical-lattice experiments. 


Entanglement is a central theoretical concept underly¬ 
ing the characterisation of quantum many-body states 
in condensed-matter and high-energy physics, as well 
as quantum information. For example, entanglement 
properties reveal exotic states of matter such as topo¬ 
logical spin liquids [T] or many-body localization mis], 
the holographic entanglement entropy identifies confine¬ 
ment /deconfinement transitions in gauge theories [H |5] , 
and entanglement is considered the central resource for 
quantum-enhanced metrology mn\ as well as quantum 
computation iHii]. In experiments, entanglement be¬ 
comes measurable via a tomographic determination of 
the many-particle quantum state (MIS], and protocols 
have been developed [16] and implemented in remarkable 
experiments m to measure entanglement entropies in 
quench dynamics and quantum phase transitions. How¬ 
ever, the resources required by these protocols scale ex¬ 
ponentially with the system size, and these experimental 
efforts are thus limited a priori to few-particle systems. 

To address the problem of detecting and quantifying 
multipartite entanglement for large systems, we consider 
below the quantum Fisher information (QFI) as an en¬ 
tanglement witness [IBHSni- Our key result is that— 
for a many-body system at thermal equilibrium at any 
temperature—the QFI can be determined directly from 
a measurement of Kubo linear response functions, in par¬ 
ticular the dynamic susceptibility (see Fig. [^. We em¬ 
phasise that this measurement prescription is indepen¬ 
dent of microscopic details of the system of interest and 
that the measurement of linear response is a standard 
tool in experiments. Importantly, only modest measure¬ 
ment resources are required that do not scale with system 
size. The presented prescription therefore makes mul¬ 
tipartite entanglement observable for a large variety of 
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Figure 1. (color online) Measurement prescription for the 
quantum Fisher information, a witness for multipartite en¬ 
tanglement. (a) A quantum many-body system is prepared 
in a thermal state at temperature T. The example shows an 
Ising spin chain (top), Eq. Q, as realisable in optical-lattice 
experiments (bottom) [21], but the general concept applies 
also to fermions and bosons, and in any spatial dimension, 
(b) Standard tools, such as inelastic Bragg or neutron scatter¬ 
ing, measure the imaginary part of the dynamic susceptibility 
x"{(^,T) as a function of frequency exchange u (green line). 
The integral Fq{T) = ^ dcjtanh (^) x^^(ct;,T) gives the 
quantum Fisher information (shaded areas), (c) This pro¬ 
cedure allows mapping out the quantum Fisher information 
as a function of temperature and transverse field 0. Low- 
temperature states (circle) can host entanglement, but it is 
lost at larger temperatures (cross). The quantum Ising chain 
has a divergence of entanglement close to the quantum crit¬ 
ical point, as well as a robust entangled region extending to 
finite temperatures (enclosed by the dashed line). 


experimental platforms, including quantum degenerate 
atomic gases as well as condensed-matter systems. 

The discovered connection between QFI and response 
functions has profound implications for our understand¬ 
ing of quantum many-body systems. For example, the 
quantum fluctuations and quantum correlations con¬ 
tained in response functions provide a rigorous lower 
bound to multipartite entanglement in thermal ensem¬ 
bles. Even more, this intimate connection to response 
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functions immediately implies universality of the QFI 
close to quantum phase transitions (QPTs), as has been 
observed theoretically in ground-state phase diagrams 
[221425] and is known for related metrics [26]. Below, 
we discuss the universal scaling laws of the QFI at van¬ 
ishing as well as non-zero temperature, and use them to 
identify a class of strongly-entangled QPTs with a diver¬ 
gent multipartiteness of entanglement and an associated 
divergent entanglement length scale EH- 

This article is organised as follows. After outlining 
background on the QFI, we present the main result of 
this work, the relation of the QFI to Kubo response 
functions. Then, we discuss its universal scaling theory, 
and discuss why QPTs yield strong divergencies while 
thermal phase transitions do not present particular 
signatures. We illustrate these completely general 
concepts by paradigmatic examples. Finally, we discuss 
some practical experimental considerations, and provide 
our conclusions. 

Background on the quantum Fisher information. — In 
recent years, the QFI has generated a lot of attention, 
because it provides a rigorous lower bound for genuinely 
multipartite entanglement [THuin]. Originally, it was in¬ 
troduced to quantify the maximal precision with which 
a parameter (a phase) 'd can be estimated using a given 
quantum state p [71 [28|. For M independent measure¬ 
ments, the quantum Fisher information Fq bounds the 
variance of 'd by (A'd)^ > l/(MFg), the so-called quan¬ 
tum Cramer-Rao bound [28] . 

Importantly, the maximal precision achievable in quan¬ 
tum phase estimation can break classical limits if the em¬ 
ployed system is in a strongly entangled A-particle state 
[7]. More precisely (see Supplementary Material [30]), if 
a state achieves a QFI density 


(where the sum includes only terms with px + Px' > 0) • 
Recently, in a remarkable atomic-gas experiment m, 
it has been demonstrated that a lower bound on the 
QFI can be measured by studying the behaviour of an 
observable’s probability distribution under the unitary 
transformation . This could then be used to 

demonstrate the presence of bipartite entanglement. 
Here, we show how the QFI can be measured directly 
and efficiently in a generic quantum many-body system 
in a thermal state at any temperature. 


Main result. — As the major result of this work, we 
rigorously relate the QFI to a Kubo response function, 

FQiT) = ^j^ dwtanh(^)x"(w,T), (4) 

where x"(ci;,T) = ^(x(co’,T)) is the imaginary, dissipa¬ 
tive part of the dynamic susceptibility in the state p with 
respect to O —the same thermal state and generator for 
which the QFI is evaluated [29] . 

Proof: The proof of Eq. 0 is straightforward. It re¬ 
quires only the minimal assumption of thermal equilib¬ 
rium. Setting h = 1 = the dynamic susceptibility is 
defined as 


X{(^,T)=iJ^ dte d{t),d ) , 


(5) 


where 0{t) = e*^^Oe“*^^. It is convenient to work in 
the Lehmann representation, i.e., the energy eigenbasis, 
where 


= F,iP\-py) 

A,A' 


(AiaiA') 
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fQ = FQ/N>m, ( 1 ) 

with m a divisor of A, then it must contain m + 1-partite 
entanglement [Hms]. 

In quantum phase estimation, the QFI determines the 
sensitivity of the state p towards a unitary transforma¬ 
tion generated by the hermitian operator O associated to 
'd. That is, it quantifies the distinguishability of p from 
p' = For a pure quantum state p = |'0)('0|, 

such as the ground state of a given Hamiltonian, the QFI 
assumes the simple form of a connected correlation func¬ 
tion, which can be easily computed or measured, 

Fq = (2) 


Matters become much more complicated in a mixed 
state, such as a thermal ensemble p = pa | A) (A|, where 
I A) is the energy eigenbasis with occupation probabilities 
Px = exp {—ExlT)jZ^ and Z is the partition function. In 
such a case, the QFI takes the considerably more complex 
structure 


A,A' 


(Pa -Pa')^ 

Pa + Pa' 


2 


(A|0|A') 


(3) 


Exploiting that for a thermal state 
2 /y dw tanh (^) 5(w + Ex' - Ex) = tanh (^^^^) = 
and correcting for prefactors, we directly obtain 
Eq. Q. This proof can be straightforwardly extended 
to the QEI matrix, and, via the fluctuation-dissipation 
theorem x"{^) = tdinh{uj/2T)S{uj), to the dynamic 
structure factor S{uj). Equation @ also presents some 
direct corollaries, such as a sum rule which we discuss in 
the Supplementary Material [30]. □ 

The conceptual importance of the identification Q is 
huge. 

Eirst of all, it makes a witness for multipartite entan¬ 
glement, the QEI, a straightforwardly measurable quan¬ 
tity. Dynamic susceptibilities are routinely measured 
in many-body systems using well-established techniques 
such as Bragg spectroscopy [ST] |32] or neutron scattering 

m- 

Second, the central Eq. 0 has fundamental theoretical 
implications. Eor example, it establishes a direct relation 
between quantum correlations contained in x"(cc;,T) to 
many-body entanglement. As Eq. shows, the entan¬ 
glement contained in x"(c(;,T) and extracted by the QEI 
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Figure 2. (color online) Universal scaling of the quantum 
Fisher information density, fq = Fq/N, calculated for the 
order parameter in the quantum Ising chain. Data at the 
critical point and for A/" = 8 ... 128. At low temperatures, the 
universal scaling laws produce an excellent data collapse. In 
a range of intermediate temperatures, the scaling approaches 
analytic predictions for the thermal critical regime in the ther¬ 
modynamic limit [34] (solid line). The strong divergence at 
small temperatures of fq ~ (dashed line) implies a di¬ 

verging multipartiteness of entanglement. 

is dominated by the thermally unaccesible high-frequency 
response. 

Third, the connection 0 has especially profound 
consequences near continuous QPTs when choosing 
for O a relevant operator in the renormalisation-group 
sense, such as the order parameter. Then, known 
universal scaling laws for x"(co’,T) translate directly 
into universal scaling for Fq. Universal scaling of the 
QFI has already been theoretically observed in ground 
states of many-body systems [221425] . In the following, 
we extend the analysis to the experimentally relevant 
regime of non-zero temperatures. 

Universal scaling of multipartite entanglement. — Con¬ 
sider a local generator O = ^ d-dimensional 

system with linear size L and N = sites. As explained 
in the Supplementary Material [30] , fq obeys the univer¬ 
sal behaviour 

fQ{T, L-\h) = L-U, . (7) 

Here, ^ is the dynamical and v the correlation-length 
critical exponent. A is the cutoff scale determined by the 
relevant perturbations 1/“^, T, and the distance from the 
critical point h. The scaling dimension of fq is Ag = 
d — 2Aq,, with Aq, the scaling dimension of Oi. 

Since fq bounds the number of entangled particles via 
Eq. its scaling behaviour allows us to identify a class 
of strongly entangled QPTs, i.e., QPTs with a divergent 
multipartiteness of entanglement. These are those tran¬ 
sitions with Ag > 0. The scaling behaviour additionally 
implies a universal length scale over which multipartite 
entanglement exists, /ent ^ fq^ ^ ^i- 2 A«/d Su.p- 
plementary Material m)- 

The framework described up to now is completely gen¬ 
eral. To illustrate its viability for understanding entan¬ 
glement in quantum many-body systems, we now turn to 
specific examples. We focus on a paradigmatic class of 


model systems presenting a QPT, namely Ising models in 
a transverse field, which are realisable in quantum-optical 
[2Q[ \2T\ as well as solid-state systems [34l |35] , 

AT N 

— = - cos 0 ^ + sin 0 ^ af . (8) 

ij=i 1=1 

Here, a = are the Pauli matrices on lattice 

sites I = 1,..., A. Depending on the Ising interactions 
Azj, this Hamiltonian exhibits a quantum critical point 
at some critical field strength Oc- The order parame¬ 
ter for the transition is {af) /N. Its strong critical 

fluctuations make the generator O = an ideal 

candidate for testing the scaling behaviour of the QFI. 

We first focus on the simplest case, the one¬ 
dimensional nearest-neighbour Ising chain, Jij = 
where the scaling exponents are known analytically, ^ = 
1 and Aq, = 1/8 [3l]. This gives Ag = 3/4, i.e., the 
nearest-neighbour quantum Ising chain lies in the class 
of strongly entangled phase transitions with divergent 
multipartite entanglement. Indeed, fq for the order pa¬ 
rameter displays a strong peak around the critical point 
Oc = 7r/4, see Fig. [^c). The entanglement radiates out 
from the peak, generating a broad entangled region also 
at non-zero T ®|37]. 

To illustrate the Ansatz Q, we consider the scaling 
with system size L and temperature T, at fixed A = 0. 
For small L, the dominant cutoff scale is A ^ L = A, 
implying 

fQ -- V3/4 , (9) 

The data in Fig. [^reproduces perfectly this strong alge¬ 
braic growth, which is remarkably close to the theoretical 
maximum of /g = A. The associated multipartite entan¬ 
glement length scale is thus highly divergent, /ent ^ A3/4. 

With increasing temperature, the cutoff scale crosses 
over to A ^ and the scaling becomes fq ^ 

This prediction can be refined by exact analytical re¬ 
sults for the dynamic susceptibility at criticality m, 
limiv^oo ^)/^ = with a function 

g that is known exactly. Performing the cj-integral in 
Eq. one immediately obtains the QEI, 

fQ{T) = C{J/TF\ (10) 

with (7 ~ 0.42. In the temperature regime of validity, 
NT IJ ^ 1 and T/J ^ 1, the exact data for finite 
chains is consistent with this scaling prediction valid 
in the thermodynamic limit (see Eig. [^. Eor T ^ J 
the system crosses over into a generic high-temperature 
asymptotic behaviour fq ^ T“^. 

Absence of signature at thermal phase transitions .— 
Remarkably, such scaling behaviour is only observed at 
quantum, but not thermal phase transitions, because 
Eq. 0 considers only quantum fluctuations. A simple 
example to demonstrate the insensitivity towards ther¬ 
mal phase transitions is provided by the fully-connected 
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Figure 3. (color online) Absence of features of the quantum 
Fisher information at thermal phase transitions, exemplified 
by the infinitely-connected Ising model, (a) Finite entangle¬ 
ment is witnessed at non-zero temperatures (enclosed by the 
dashed line), and a divergence appears near the quantum crit¬ 
ical point at 0 = 7r/4 and T = 0. But, crucially, /q shows 
no features at the thermal phase transition (solid line). Data 
for N = 1000. (b) Derivatives of fq remain also smooth at 
the thermal transition (vertical line). Top: fq (dashed) and 
dfq/d{T/J) (solid line). Bottom, from dark to light: second 
to fifth derivative. Data for N = 1600 and 0 = tt/S. 


transverse-field Ising model, Jij = l/N, V/, j, similar to 
the model describing the experiments of Ref. [20]. In 
contrast to its nearest-neighbour counterpart, this model 
exhibits, additionally to the QPT at Oc = 7r/4, also a 
thermal phase transition [38] . 

Figure |^a) shows fq for the order parameter in the 
temperature-transverse-field plane. We delegate its scal¬ 
ing analysis to the Supplementary Material m- More 
important at this point, while fq shows a divergence at 
the QPT, no particular feature can be discerned at the 
thermal phase transition. Neither do such features ap¬ 
pear in derivatives of fq [Figure |^b)]. 

However, it is known that the static isothermal 
susceptibility of Ihe order parameter diverges [39] . 
The reason why the QFI instead remains featureless 
becomes clear when decomposing = Xei + XvV iiilo 
its two fundamental parts m, the elastic (or Curie) 
contribution Xeb and the quantum-mechanical van-Vleck 
correction XvV, which is continuously connected to the 
Kubo susceptibility, XvV = ^)- H is Xei 

that diverges at a thermal phase transition. As is shown 
in the Supplementary Material m, Xei can be related 
to the Fisher information in a classical scenario [41] that 
has no relation to entanglement. The second term XvV^ 
on the other hand, remains smooth at thermal transi¬ 
tions. The QFI, thus, considers only the contribution to 
the susceptibility that is due to quantum fluctuations 
and remains insensitive to thermal phase transitions. 

Experimental considerations .— Let us finally address 
some practical aspects that will be important for exper¬ 
iments. The measurement prescription Q is very flexi¬ 
ble, since its proof did not make any assumptions on mi¬ 
croscopic details of the system under study (other than 
thermal equilibrium). As a consequence, it applies in 
any spatial dimension, for any hermitian generator O, 


and it can be equally used for systems of spins, bosons, 
or fermions. For example, the quantum Fisher informa¬ 
tion may also be probed in ongoing atomic-gas experi¬ 
ments that aim at studying low-temperature phases of 
the Fermi-Hubbard model [42] 03]. A bound for mul¬ 
tipartite spatial entanglement as in Eq. 0 is known to 
exist in all cases where the generator is a sum of local 
operators, O = Oi, when Oi has a bounded spectrum 
[7] (see Supplementary Material [30]). 

Importantly, the complexity of measuring 
does not scale with system size, thus allowing for an 
efficient evaluation of the QFI via Eq. Q. Moreover, 
the scaling analysis overcomes a usual practical difficulty 
for studying the QEI, the optimal choice for the linear 
generator O = '^lOi. In the vicinity of a QPT, the 
choice becomes clear: one may select any suitable, rel¬ 
evant operator O/, preferably the one with the largest 
scaling exponent Ag. A good choice will often be the 
order parameter, as in the examples above. The Supple¬ 
mentary Material m contains an example for the Mott- 
insulator-superfluid transition, where a generator differ¬ 
ent from the order parameter but accessible by Bragg 
spectroscopy, allows one to extract universal behaviour. 

Eurther, the prescription is robust against intrinsic 
sources of imperfections. Eor example, an uncertainty in 
determining the system’s temperature can be mitigated 
by choosing the lowest temperature estimate. Due 
to the monotonicity of the tanh(co’/2T) as a function 
of frequency, this produces a lower bound for the 
QEI. Another natural error source is a finite spectral 
resolution. Assuming a symmetric spectral broadening, 
the concavity of the tanh assures that the integration 
again returns a lower bound to the true QEI, with less 
severe results at smaller temperatures. Also, the integral 
in Eq. 0 will be limited to some frequency range (the 
upper integral limit represents, as usual, a scale much 
larger than energies accessible by the considered degrees 
of freedom). This limitation will again deliver a lower 
estimate. Importantly, non of these errors will produce 
a false positive indicator for entanglement. Even more, 
close to a QPT, the existence of the universal scaling 
laws assures an inherent robustness of the QEI. There¬ 
fore, in contrast to some engineered highly-entangled 
states, such as the GHZ state, the entanglement wit¬ 
nessed by the QEI close to QPTs is unaffected by weak 
symmetry-preserving perturbations. 

Conclusion. — To conclude, we have developed a sim¬ 
ple prescription for measuring the quantum Eisher in¬ 
formation, which is independent of spatial dimensional¬ 
ity or other microscopic details, and relies only on well- 
developed tools for measuring dynamic susceptibilities. 
Many-particle entanglement thus becomes observable in 
various experimental platforms, including quantum de¬ 
generate atomic gases as well as solid-state samples. 

The discovered prescription adds an important tool 
to the ongoing search for measurable entanglement wit¬ 
nesses gn Ss], some of which have already allowed ex- 
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periments to demonstrate the presence of bipartite en¬ 
tanglement [20l l46f[49] . The prescription is also comple¬ 
mentary to efforts for measuring the entanglement en¬ 
tropy which can reveal more complete information, 
at the expense, however, of an exponential increase in 
measurement resources. Moreover, generic scaling be¬ 
haviour of the block entanglement entropy is rigorously 
demonstrated only near critical points that are described 
by a conformal field theory [44]. In contrast, the QFI dis¬ 
plays scaling behaviour in the vicinity of arbitrary contin¬ 
uous quantum phase transitions. This scaling behaviour 
allowed us to classify quantum phase transitions with a 
divergent multipartiteness of entanglement. 

In order to achieve an advantage over classical systems, 
identifying and quantifying large-scale entanglement 
in many-particle systems will be an indispensable 
ingredient for quantum simulation EHII] and quantum 
metrology m n\- It will therefore be worthwhile for 
future studies to generalise the concept developed here 
to other kinds of phase transitions, e.g., topological, or 
to states beyond thermal ensembles, such as generalised 
Gibbs ensembles. 


Supplementary material 

In this supplementary material, we (i) provide addi¬ 
tional technical details, (ii) explain in detail the scal¬ 
ing behaviour for the quantum Fisher information and 
present additional numerical results, (hi) discuss the rela¬ 
tionship Fisher-information/susceptibility in a classical, 
thermal setting, and (iv) present a sum rule that may be 
used to bound the quantum Fisher information. 


SI. RELATION OF QFI TO MULTIPARTITE 
ENTANGLEMENT 

The QFI has a deep connection with multipartite en¬ 
tanglement [iHiiin]- Consider a system of N spins with 
length S' = i, and a linear observable O = ^ X]/ , 

where ai = (erf , erf , af) is a vector containing the Pauli 
matrices erf associated to spin /, and = (nf,nf,nf) is 
a unit vector on the Bloch sphere. Then the system hosts 
at least m + 1-partite entanglement if the QFI associated 
to O fulfils [18] 
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Fq > 


N 

m 


m 


N ■ 


N 

m 


m 


(SI) 


where [X\ is the largest integer smaller than or equal 
to X. Eor m a divisor of A, the condition ( |S1| ) attains 
the simple form Q when expressed through the quan¬ 
tum Eisher information density /g = Fq/N. Note that 
typically response functions that measure linear opera¬ 
tors contain a factor of 1/Y relative to our definition, 
Eq. (§, so that they will observe directly /g rather than 

Fq. 

The proof of Ref. [18] for bounding the multipartite 
entanglement with Eq. (SI) can be directly translated 


to degrees of freedom other than spins 1/2, as long as 
O represents a sum of local operators with bounded 
spectrum [7]. If hmax and hmin denote the largest, 
respectively smallest, eigenvalue of O, then the right 
hand side of the condition (SI) acquires the prefactor 
(^max — ^min)^- Therefore, the QEI can also witness 
spatial entanglement in systems other than spins-1/2, 
such as larger spins or fermions. Note that the relation of 
the QEI to response functions, Eq. is independent of 
any such microscopic details of the underlying quantum 
many-body system or the hermitian operator O, which 
may even be non-local or unbounded. 


S2. SOLVABILITY OF THE CONSIDERED 
MODELS 

All models used in this article for illustrating the main 
concepts are exactly solvable. The one-dimensional Ising 
chain in a transverse field [Eq. ^ with Xij = 
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can be mapped to a free-fermion problem [34]. Dynami¬ 
cal susceptibilities can then be calculated via Wick’s de¬ 
composition of expectation values m, and the Pfaffians 
appearing in the resulting expressions can be evaluated 
efficiently using the algorithm described in Ref. m- 
The infinite-range Ising Hamiltonian [Eq. Q with 
Jij = 1/N] commutes with both + 

and where erf/2, a = z. As a con¬ 

sequence, the Hamiltonian decomposes into disconnected 
blocks when represented in the common eigenbasis of 
and . Each block grows linearly with particle number 
N and can be diagonalized efficiently. Here, we consider 
the largest of these blocks, with dimension N 1. The 
infinite-range Ising model has a thermal phase transition 
for ^ < 7r/4, with critical temperature [38] 


J log[(l + tan(6>))/(1 — tan(6>))] 


at low energy, i.e., they display the same power-law de¬ 
cay of correlations, which is described by a small set of 
universal exponents { }. 

In finite systems, or systems at finite temperature, the 
power-law of the correlations is cut off by system size or 
the thermal correlation length. In such cases, one can 
extract the universal properties of the QPT by perform¬ 
ing an appropriate analysis of scale transformations, i.e., 
transformations that modify the lattice spacing a ^ Aa, 
where A > 1 (this scale transformation can be envisioned 
as grouping subsets of A spins together). Close enough to 
the critical point one can expand all local operators in the 
basis of scaling operators, i.e., those that transform un¬ 
der rescaling with a simple power of A, 
with scaling exponent Aq, > 0. The two-point correlation 
function of these operators decays as 

(C>“(0)C>“(r)) - (C»“(0))(d“(r)) oc . (S3) 


In all models, we set /i = = a = 1, where a is the 

lattice spacing. 


S3. SCALING AND UNIVERSALITY 


In the main text, we have demonstrated a close re¬ 
lationship of the QEI to thermodynamic response func¬ 
tions. These response functions underlie universal scal¬ 
ing behaviour close to quantum critical points, and con¬ 
sequently directly imprint a universal scaling behaviour 
onto the QEI. Aspects of critical scaling of the QEI and 
related quantum metrics have also been discussed, e.g., 
in Refs. [22ll26] . The purpose of the present supplemen¬ 
tary section is to present a general, unified scaling theory 
for the QEI. 


A. Basic principles 

Before discussing the scaling of the QEI in detail, let 
us briefly recall the basic principles of the scaling Ansatz 
m- Here, we focus on lattice models with lattice spacing 
a that have a continuous phase transition at zero tem¬ 
perature, viz a quantum phase transition (QPT). Gen¬ 
eralisation to continuum models is straightforward. The 
dimensionality of the system is d, its linear size L, and 
the total particle number N = (where we set a = 1). 
To simplify the discussion, we focus on the simplest sce¬ 
nario, QPTs obtained by tuning a single relevant field, 
say the magnetic field h, to the critical point he- In this 
scenario, any small variation of h = \h—hc\ drives the sys¬ 
tem in the thermodynamic limit to its stable fixed point 
with zero correlation length. Nevertheless, in the scaling 
regime, i.e., sufficiently close to the phase transition, all 
low-energy physics becomes universal. As a remarkable 
consequence of this universality, very diverse microscopic 
models develop the same collective emerging behaviour 


In particular, the fixed-point Hamiltonian is a scaling 
operator with scaling exponent 

Scaling operators that transform with Ac^ > /\h are 
called irrelevant operators, since, when added as a weak 
perturbation, coarse graining decreases their importance 
relative to the original Hamiltonian. Those operators 
with Ack < Aj:/, on the other hand, are called relevant^ 
since they become larger and larger under coarse-graining 
transformations, eventually driving the system to a dif¬ 
ferent fixed point. In the sealing regime^ the regime where 
relevant perturbations are still small compared to the 
original Hamiltonian, one can extract the universal low- 
energy physics by studying the response of the system to 
scaling transformations. 

Both the inverse system size L~^ and the finite temper¬ 
ature T are relevant perturbations. They increase under 
a rescaling as L~^ \L~^ and T A^T, where 2 : is 

the dynamical critical exponent. The scaling regime for 
a finite and cold system is thus defined by <C 1 and 
T ^ 1. We will now discuss the scaling behaviour of the 
QEI with respect to these perturbations, as well as to a 
deviation h from criticality. 


B. Scaling of the quantum Fisher information 


To construct a scaling theory for the QEI close to a 
quantum critical point, we start from Eq. Q, which de¬ 
fines the QEI at zero temperature, and consider a local 
generator O = Of . By inserting Eq. ( [S^ into Eq. Q , 
one finds that /g transforms under rescaling by a factor 
A as 


/q(T', L’-\h') = L-^A, + C. 

(S4) 

Here, C is a constant that is related to the non-universal 
short-distance correlations, and n is the critical expo¬ 
nent associated to the correlation length. According 
to Eq. (S4), the QEI has the scaling exponent Ag = 
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(i —2 Aq,. The scaling function (/)q(TA^, en¬ 

codes the large-distance part of the correlation functions, 
which dominates the scaling if the condition Ag > 0 is 
fulfilled. In the following discussion, we assume that this 
is the case, allowing us to neglect the non-universal con¬ 
stant C, leading to Eq. 0. Below we will encounter, in 
the guise of the hard-core boson chain, an example where 
Aq ^0- In such cases, derivatives of /g still allow one 
to extract the purely universal contribution. Note that 
when Ag = 0 there can be logarithmic corrections to 
Eq. (ISl. 


We can use the scaling form Eq. (S4) as long as all 
relevant perturbations are small, that is as long as we 
stay in the scaling regime. The most relevant perturba¬ 
tion will constitute the most significant breaking of scale 
invariance, a fact that enables us to extract valuable in¬ 
formation about the involved critical exponents. Eor ex¬ 
ample, when both L~^ and T are sufficiently small, the 
deviation h from the critical point represents the most 
relevant perturbation: Scale invariance will be broken at 
the scale A ^ h~^. This implies that we can trade A by 
and varying h in this regime has the same effect as 
a rescaling transformation. In this regime, we expect /g 
therefore to behave as 


/q(T, L-\h) = {Th--\ L-^h-^) . (S5) 

Eor sufficiently relevant operators, for which Ag > 0, 
/g thus diverges as when approaching the criti¬ 

cal point. In the above expression, (j)Qi(Th~^^^ 
is a scaling function, embodying the fact that the nor¬ 
malised quantum Eisher information density, ^ 

is a well defined function of the scaling variables Th~^^ 
and rather than of L, T, and h separately. 

Upon decreasing h, the system will reach a crossover 
with the next relevant perturbation. If L-i < the 

crossover will occur around ~ h'' and carry over 

into a thermally dominated regime, for which the scaling 
of the QEI is dictated by 


(S6) 

This scaling behaviour characterises the robustness of 
critical entanglement against non-zero temperatures (see, 
e.g., Ref. [S] for a discussion of the robustness in presence 
of noise). 

In the opposite case L -1 > 

upon reducing h one 
enters a regime dominated by finite-size effects. Here, /g 
behaves as 


fQ(T,L-\h) = (TL^L^/’'h}j . (S7) 


In this regime, /g saturates to a finite, L dependent 
value, which diverges with system size as . The 
crossover is located around L~^ ^ and should thus 
approach the critical value h = 0 as . Another in¬ 

teresting scenario appears when T and 1/L are compara¬ 
ble, giving rise to further complex crossover phenomena. 


This can be seen in Eig. of the main text, where we 
study the universal scaling behaviour of the QEI in the 
paradigmatic model of the quantum Ising chain. 

This scaling analysis immediately permits us to iden¬ 
tify strongly entangled quantum critical points. These 
are those critical points where Ag > 0, i.e., where at low 
temperatures and close to criticality, /g diverges with 
system size. Additionally, the scaling also allows us to 
extract a length scale of the entanglement present in the 
system EH- Combining Eq. ( |S4| and main-text inequal¬ 
ity 0, one sees that the number of entangled particles 
scales as Nq > /g ^ ^d-2Ao, ^ which defines a multipar¬ 
tite entanglement length scale 


lent ~ . (S8) 

Here, we assumed an isotropic distribution of entan¬ 
glement; for anisotropic situations, the longest distance 
along which particles are witnessed to be entangled is ac¬ 
tually larger. This entanglement length is related to the 
correlation length ^ of the generator O. Since at crit¬ 
icality the correlation length scales as ^ ^ A, one has 
^ent ^ ^i- 2 AQ,/d^ 2Aa/d < 1, the length scale over 
which particles are entangled with each other thus di¬ 
verges with system size. 


C. Universal scaling of in the 
Mott-insulator—superfluid transition 


The scaling behaviour of /g, Eq. (S4), appears in all 
relevant operators, not only the order parameter. As 
an illustration for this fact, we choose a system that is 
relevant for experiments on cold gases, a system of ultra¬ 
cold bosonic atoms confined in a one-dimensional optical 
lattice. In the limit of strong on-site interactions, this 
system can be modelled by the Hamiltonian [52] 


-^hcb — 


AT-l N 

E , (S9) 


1=1 


1=1 


where the bosonic creation (annihilation) operators bj 
(bi) obey the hard-core constraint blbj = 0. This model 
has a QPT at fic/J = 1 from a critical superfluid phase 
to a Mott-insulator at uniform filling, with critical expo¬ 
nents z = 2 and u = 1/2 [34] . 

The transition has several peculiarities. Eirst, the ex¬ 
act ground state is fully separable at the QPT [53]. Sec¬ 
ond, one of its order parameters, the particle density 
h = ^ibjbi/N^ commutes with the Hamiltonian. We 
can, however, use a modified observable that is still a rel¬ 


evant operator, the staggered density O = ^li—'^ybjbi. 
In optical-lattice experiments, the dynamic susceptibil¬ 
ity of such density fluctuations may be probed by Bragg 
spectroscopy [nisi]. 

The Hamiltonian Eq. (89) is exactly solvable with the 
same methods as for the Ising chain in a transverse field 
ISl El. Moreover, under periodic boundary conditions. 






(a) 2 
TU 


4n /i < 0 

4(1 — n) /i > 0 



Figure SI. (color online) Quantum Fisher information den¬ 
sity, /q, for the staggered density of hard-core bosons in a 
one-dimensional optical lattice, (a) In part of the super¬ 
fluid phase, fq witnesses two-site entanglement. Data for 
N = 1000. (b) Inset: dfqjdii shows a peak that, with de¬ 
creasing T, diverges and moves towards the critical point. 
Main panel: The predicted scaling exponents produce a per¬ 
fect data collapse for dfq/dfx. Data for N — 100 and, from 
blue to red, 8 different temperatures distributed on a loga¬ 
rithmic scale between T = 0.01 — 0.2J (inset: 4 temperatures 
in 0.01 -0.13J). 


the dynamic structure factor S(^q^ lS) associated to density 
fluctuations O = bj h at wave vector q can be 

evaluated analytically [55]. It is gapless at g = 0 and 
acquires a maximum at g = tt, where 

S{'K,uj) = 27r^(5(w -ek + efe+^)/(efe+^)[l - /(e^)], 

(SIO) 

with e/e = —2[J + /icos(/c)] the single-particle en¬ 
ergy of the diagonalised Hamiltonian and f{^k) = 

l/[exp(e/e/T) + 1] the Fermi-Dirac occupation. As a 
consequence of the gaplessness, the QFI evaluated at 
g = 0 vanishes, while it assumes its maximum value 
at g = TT, i.e., for the staggered density. Using the 
fluctuation-dissipation theorem and evaluating the inte¬ 
gral in Eq. the corresponding quantum Fisher infor¬ 
mation density becomes 

= V f{ek+-^)[l - f{€k)]. 

k ^ ' 

^ (sii) 

This quantity is plotted in Fig. [Sl[a) as a function 
of temperature and chemical potential, for N = 1000. 
The entanglement witnessed by the QFI extends up to 
T ^ 0.91J, improving over other witnesses showing en¬ 
tanglement up to T « 0.7J [511 EH- In contrast to 
Ref. [56], however, the QFI for the staggered magneti¬ 
sation does not witness genuine multipartiteness of en¬ 
tanglement, but saturates at the value /q = 2. (Notably, 
the QFI for other generators such as hi + h\ does show 
divergent behaviour.) Towards the phase transition, the 
QFI for the staggered density decreases, consistent with 
the presence of the factorisation point. Despite this ab¬ 
sence of divergences, the phase transition imprints its 
character onto the universal behaviour of the QFI and 
its derivatives, as explained in what follows. 

At T = 0 and N = oc, one can use Eq. <© to derive a 
particularly simple analytical formula for the QEI of the 


staggered density. 


(S12) 


Here, we defined the mean particle density n = (n), 
which is n = 1 — acos(/i/J)/7r for |/i| < J and 0 else. 
The quantum Eisher information density attains its max¬ 
imal value of /q = 2 at /i = 0 and decreases towards 
the phase-transition point, where it vanishes. The above 
analytic expression allows us to obtain the scaling di¬ 
mension of fq by expanding it close to the critical point, 
fq oc with h = 1 — /i/J. Comparison with the 

scaling form Eq. (S5) gives the exact scaling exponent 


Aq = — 1. Since Aq < 0, /q of the staggered magneti¬ 
sation displays no divergent behaviour. 

Nevertheless, as seen in the inset of Eig. [STj^b), the 
non-analyticity of /g at = 0 leads to a divergence in 
its derivatives. Considering the temperature as the most 
relevant perturbation, Eq. (S6), we obtain the universal 
behaviour 


_ j^-Aq I Z + ll {ZV) 

dh 


(S13) 


That is, dfq/dh has a peak that diverges as 
j^-^qIz+iI{zv) position of which lies at a line 

h ^ Using the analytically known exponents 

—Aq/z^l/{zTy) = 1/2 and l/{zu) = 1, an excellent scal¬ 
ing collapse is achieved up to rather large temperatures 
[main panel of Eig. |Sl[b)]. We also performed fits to the 
numerical data for the peak height and peak position. 
These yield exponents of 0.49 and 1.02, respectively— 
even without high-precision data, we are able to repro¬ 
duce the exact scaling exponents within two percent ac¬ 
curacy. 


D. Universal scaling in the infinite-range quantum 
Ising model 

In the main text, we considered the infinite-range 
Ising model in a transverse field, main-text Eq. ^ with 
Jij = 1/A, to illustrate the absence of features in the 
QFI at thermal transitions (see Fig.[^. Here, we provide 
its scaling analysis in the vicinity of its quantum phase 
transition. 

The infinite-range Ising model describes the universal 
properties above the upper critical dimension, with the 
longitudinal magnetization O = order 

parameter. In Fig. we plot the QFI as a function 
of temperature for different system sizes at the system’s 
quantum critical point, = 0. Using that Ag = 1/3 
and 2 : = 1/3 [58] one obtains the scaling of the quantum 
Fisher information density fq ^ A^/^0(TA^A), From 
a scaling collapse, shown in Fig. we determined nu¬ 
merically that Ag =0.34 and 2 : = 0.33 which is in very 
good agreement with the exact solution. 
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Figure S2. (color online) Universal scaling of the quantum 
Fisher information density, fq = Fq/N, calculated for the 
order parameter in the quantum Ising model with infinite- 
range interactions. Finite-temperature and finite-size scaling 
at criticality 0 = 7r/4. Bullets are from lighter to darker blue 
for iV = 100,, 3000. The universal scaling laws produce an 
excellent data collapse. The divergence fq r\J jVl/3 

is much 

weaker than in the one-dimensional quantum Ising model, 
reflecting the more mean-field-like character of the infinite- 
range model. Nevertheless, it still implies a diverging multi¬ 
partiteness of entanglement. 


S4. CLASSICAL CONTRIBUTION TO QFI AND 
BEHAVIOUR AT THERMAL PHASE 
TRANSITIONS 


In the main text, we study how well thermal states that 
lie close in a phase diagram are distinguishable, focusing 
entirely on quantum effects. More precisely, we ask how 
fast the probability distribution of a thermal state p = 
^xPx 1^) (^1? with Pa = exp(—F^ a/T)/Z, changes due to 
the unitary transformation with . It is instructive 
to extend this analysis to any source of distinguishability, 
be it quantum or classical. 

Consider a canonical ensemble p at temperature T for 
a Hamiltonian H. We wish to compare this state with the 
equilibrium ensemble p' at the same temperature but for 
H' = where 6P 1 and where J represents a 

suitable energy scale. There are two sources that render 
p and p' distinguishable. The first is the one considered 
in the main text. It can be rephrased as an adiabatic 
evolution of the state by slowly changing the Hamilto¬ 
nian from H to H'. If the evolution is adiabatic and 
unitary, the populations px remain unchanged. Any dif¬ 
ference between initial and final density operator is then 
solely due to a modification of the energy eigenstates. 
Taking a unit time step t = 1/J^, the change amounts to 
I A) ^ 1 ^^^ This is formally equivalent to the phase 

estimation scenario considered in this article. 

However, one may extend the scenario by including a 
sufficiently long waiting time after the transfer, such that 
the populations at the new Hamiltonian H' can rether- 
malise with an external bath. This change of populations 
provides a second source of distinguishability between p 
and p', which to leading order remains disparate from 
the one described above. We can incorporate this effect 


by considering the full QFI [28] 


Fq = Fc F Fq 

^ y. {dpx/ddf ^ 2 V -P\'? 


(S14) 

(A|0|A') ' . 


The second term, Fq, the subject of this work, is the 
term that is relevant at a QPT at T = 0, where the first 
term, vanishes identically. 

Both terms are related to different aspects of the sus¬ 
ceptibility of p towards the generator O. As demon¬ 
strated in Eq. (2) of the main text, Fq is in one-to-one 
correspondence to the Kubo dynamic susceptibility at 
non-zero frequencies. In Kubo’s linear response theory, 
a system is perturbed weakly in a time-dependent way, 
generating a unitary dynamics, but—to leading order— 
without changing the populations of energy levels m- 

The change of probabilities, instead, is related to the 
static, isothermal susceptibility, as has already been 
pointed out within purely classical statistical physics 
m- Here, we demonstrate this relation in a quantum 
statistical calculation. Consider a static perturbation, 
H ^ H' = H F SP'O, with = T6P (i.e., we choose 
the energy scale J = T). The expectation value of the 
observable O is given by a derivative of the free entropy 
InZ, with Z = ^^exp(—Fy/T) the partition function. 


(O) =T 


d\nZ 

dd' 




dEx 

d'd' ’ 


(S15) 


and the corresponding static, isothermal susceptibility is 

(S16) 


^ (^) _ ^Px dEx d^Ex 


The derivative of px reads 

dpx Px / dEx 




dP' 


dP' 


-E 


dEy 

d'd' 


(S17) 


One can use perturbation theory to evaluate the change 
of energy levels following the static perturbation. Work¬ 
ing in a basis where degenerate eigenlevels |A) are predi- 
agonalised with respect to the perturbation operator (9, 
we obtain 


= Xel + XvV • 



\dx,x'?{px-p'x) 

Ex' - Ex 


(S18) 


Here, we defined the abbreviation Ox,x' = The 

first term is the elastic contribution, or classical Curie 
contribution, arising from a change of populations. It 
vanishes in a pure state, such as the quantum ground 
state at T = 0. The second term is the quantum- 
mechanical van-Vleck correction, arising from matrix el¬ 
ements between different energy levels. It is connected 
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smoothly to the real part of the Kubo dynamic suscep¬ 
tibility at non-zero frequencies and vanishes in a purely 
classical theory. Importantly, the elastic peak is a sin¬ 
gular contribution at exactly cj = 0, and does thus not 
influence Pq- Since it is this elastic peak that diverges 
at thermal phase transitions, these do not appear as par¬ 
ticular features in Fq- 

We can connect the elastic peak to the classical con¬ 
tribution to the Fisher information [41]. Equation (S17) 

yields Fq = Ea = Xei/T. Similarly to Fq 

thus has two parts: Fq, the subject of this article, re¬ 
lated to entanglement and divergent at QPTs; and Fc, 
of non-unitary, classical origin and divergent at thermal 
phase transitions. Both are related to the order param¬ 
eter’s susceptibility. Remarkably, by choosing the sus¬ 
ceptibility at vanishing or non-vanishing frequencies, we 
can extract the purely classical or the purely quantum- 
mechanical contribution to the quantum Fisher informa¬ 
tion. 


S5. SUM RULE 

Our central relation Eq. (2) allows us to derive a 
Thomas-Reiche-Kuhn-type [59l[6T] sum rule that may be 
useful to bound the QFI. Since tanh (^) < we have 

1 2 7 ^ 

Fq{T)<-- / dc^c^x"(c.,T). (S19) 

-I- ^ Jo 


Due to the Kramers-Kronig relations, we can now take 
advantage of a superconvergence theorem [62]. This 
theorem states that two complex functions g and /, 
for which g{x) = P dx holds, are related via 

dxf{x) = lim^^^oo {x‘^g{x)) , with P the principal 
part of the integral. By setting f{x) = —^xx"{x) and 
g{x) = x'{x) = 5?(x(x)), one obtains 


Fq{T) < lim {u\'{u,T)) , (S20) 

J_ uj—yoo 


a Thomas-Reiche-Kuhn-type sum rule that upper 
bounds the QFI. Here, as in Eq. of the main text, the 
limit to infinity is to be understood as a frequency range 
much larger than all energy scales of the considered 
degrees of freedom (but smaller than energy scales 
at which resonances to additional degrees of freedom 
appear, the entanglement between which, however, could 
be analysed analogously). It is interesting to note that, 
if we t ake t he classical limit by sending T ^ oo, the 
bound ds^ becomes an equality relating the QFI to 
the high-frequency limit of x'i^)- 
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